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Abstract 

We briefly review certain aspects of cosmic microwave background anisotropies as 
generated in passive and active models of structure formation. We then focus on 
cosmic strings based models and discuss their status in the light of current high- 
resolution observations from the BOOMERanG, MAXIMA and DASI collabora¬ 
tions. Upcoming megapixel experiments will have the potential to look for non- 
Gaussian features in the CMB temperature maps with unprecedented accuracy. We 
therefore devote the last part of this review to treat the non-Gaussianity of the 
microwave background and present a method for computation of the bispectrum 
from simulated string realizations. 
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1 Introduction 


Anisotropies of the Cosmic Microwave Background radiation (CMB) are di¬ 
rectly related to the origin of structure in the universe. Galaxies and clusters 
of galaxies eventually formed by gravitational instability from primordial den¬ 
sity fluctuations, and these same fluctuations left their imprint on the CMB. 
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Recent balloon ®§ and ground-based interferometer ([|) experiments have 
produced reliable estimates of the power spectrum of the CMB temperature 
anisotropies. While they helped eliminate certain candidate theories for the 
primary source of cosmic perturbations, the power spectrum data are still 
compatible with theoretical estimates of a relatively large variety of models, 
such as ACDM, quintessence models or some hybrid models including cosmic 
defects. These models, however, differ in their predictions for the statistical 
distribution of the anisotropies beyond the power spectrum. The MAP (cur¬ 
rently in space) and Planck (scheduled for launch in 2007) satellite missions 
will provide high-precision data allowing definite estimates of non-Gaussian 
signals in the CMB. It is therefore important to know precisely what are the 
predictions of all candidate models for the statistical quantities that will be 
extracted from the new data, and to identify specific signatures of the various 
models. 

There are two main classes of models of structure formation —passive and ac¬ 
tive models. In passive models, density inhomogeneities are set as initial con¬ 
ditions at some early time, and while they subsequently evolve as described by 
Einstein-Boltzmann equations, no additional perturbations are seeded. On the 
other hand, in active models a time-dependent source generates new density 
perturbations through all time. 

Most realizations of passive models are based on the idea of inflation. In sim¬ 
plest inflationary models it is assumed that there exists a weakly coupled 
scalar held 0, called the inhaton, which “drives” the (approximately) expo¬ 
nential expansion of the universe. Quantum fluctuations of cj) are stretched by 
the expansion to scales beyond the horizon, thus “freezing” their amplitude. 
Inflation is followed by a period of thermalization, or reheating, during which 
energy is transferred to usual particles. Because of the spatial variations of </> 
introduced by quantum fluctuations, thermalization occurs at slightly different 
times in different parts of the universe. Such fluctuations in the thermalization 
time give rise to density fluctuations. Because of their quantum nature and 
because of the fact that initial perturbations are assumed to be in the vac¬ 
uum state and hence well described by a Gaussian distribution, perturbations 
produced during inflation are expected to follow Gaussian statistics to a high 
degree, or either be products of Gaussian random variables. This is a fairly 
general prediction that will be tested shortly with MAP and more thoroughly 
in the future with Planck data. 

Active models of structure formation are motivated by cosmic topological de¬ 
fects with the most promising candidates being cosmic strings. It is widely 
believed that the universe underwent a series of phase transitions as it cooled 
down due to the expansion. If our ideas about grand unification are correct, 
then some cosmic defects, such as domain walls, strings, monopoles or textures, 
should have formed during phase transitions in the early universe (3j). Once 
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formed, cosmic strings could survive long enough to seed density perturba¬ 
tions. Generically, perturbations produced by active models are not expected 
to be Gaussian. 

Defect models possess the attractive feature that they have no parameter free¬ 
dom, as all the necessary information is in principle contained in the underlying 
particle physics model. This is a feature that fascinated Dennis Sciama: during 
his days at Trieste, he would always be willing to discuss these “marvelous cos¬ 
mic defects” with everyone. He even encouraged a few students sans liaisons 
contraignantes to pursue studies on this promising subject and facilitated all 
necessary means for it. The relevant, ubiquitous dimensionless quantity for 
every astrophysical and cosmological signature left by strings is G/i, where G 
is Newton’s constant and /i stands for the mass per unit length of the string, 
with a value of roughly G/r ~ 10~ 6 for GUT strings. This happens to be just 
the correct order of magnitude of the level of CMB anisotropies. Back in 1996, 
the power spectrum plot of 1(1 + 1)G/ vs. I consisted of some scattered points 
and strings were considered a likely candidate for the origin of structure for¬ 
mation. Moreover, the error bars were so big that Dennis would never miss 
the opportunity to tease inflation fans who wanted to see a peak in that plot: 
“Come on”, he would say, “you do not claim you actually see a peak there, 
do you..!”. 

To close this section, let us mention that, in addition to purely active or 
passive scenarios, perturbations could be seeded by some combination of the 
two mechanisms. For example, cosmic strings could have formed just before 
the end of inflation and partially contributed to seeding density fluctuations. 
It has been shown (^) that such hybrid models are very successful in fitting the 
CMB power spectrum data. Cosmic strings are generally expected to produce 
distinguishing non-Gaussian features in the CMB and it will soon become 
possible to look for them in the data from MAP and Planck. 


2 CMB power spectrum from strings 


Calculating CMB anisotropies sourced by topological defects is a rather dif¬ 
ficult task. In inflationary scenario the entire information about the seeds is 
contained in the initial conditions for the perturbations in the metric. In the 
case of cosmic defects, perturbations are continuously seeded over the pe¬ 
riod of time from the phase transition that had produced them until today. 
The exact determination of the resulting anisotropy requires, in principle, the 
knowledge of the energy-momentum tensor [or, if only two point functions are 
being calculated, the unequal time correlators (|j)] of the defect network and 
the products of its decay at all times. This information is simply not available! 
Instead, a number of clever simplifications, based on the expected properties 
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Fig. 1. Data from BOOMERanG (solid error bars), MAXIMA (dotted error bars) 
and DASI (dashed error bars). The line shows the CMB power spectrum produced 
by the wiggly local string model of Ref. @ in a closed universe with D tota i = 1.3, 
^baryon = 0.05, Rcdm = 0.35, = 0.9, and Hq = 65 km s^ 1 Mpc _1 . 

of the defect networks ( e.g . scaling), are used to calculate the source. The 
latest data from BOOMERanG ([[]), MAXIMA (j2|) and DASI @ experiments 
clearly disagree with the predictions of these simple models of defects (^). 

The shape of the CMB angular power spectrum is determined by three main 
factors: the geometry of the universe, coherence and causality. 


The curvature of the universe directly affects the paths of light rays coming to 
us from the surface of last scattering (||) . In a closed universe, because of the 
lensing effect induced by the positive curvature, the same physical distances 
between points on the sky would correspond to larger angular scales. As a 
result, the peak structure in the CMB angular power spectrum would shift to 
the larger angular scales or the smaller values of l. 


The prediction of the cosmic string model of (|9|) for fl tota i = 1.3 is shown in 
Fig. 1. As can be seen, the main peak in the angular power spectrum can be 
matched by choosing a reasonable value for k2 t otai- However, even with the main 
peak in the right place the agreement with the data is far from satisfactory. 
The peak is significantly wider than that in the data and there is no sign of a 
rise in power at l ~ 600—800 suggested by the data. The situation is even worse 
for global strings which fail to produce a distinct peak of any kind. For local 
strings, the sharpness and the height of the main peak in the angular spectrum 
can be enhanced by including the effects of the gravitational radiation (|T0|) 
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and wiggles (|9]). More precise high-resolution numerical simulations of string 
networks in realistic cosmologies with a large contribution from Da are needed 
to determine the exact amount of small-scale structure on the strings and the 
nature of the products of their decay. It is, however, unlikely that including 
these effects alone would result in a sufficiently narrow main peak or any 
presence of secondary peaks^J. It seems that now Dennis would not have had 
strong arguments to tease his friends. 

Let us next consider issues of causality and coherence and how the random 
nature of the string networks comes into the calculation of the anisotropy 
spectrum. 

Both experimental and theoretical results for the CMB power spectra involve 
calculations of averages. To estimate the correlations of the observed tem¬ 
perature anisotropies, experimentalists compute the average over all available 
patches on the sky. When calculating the predictions of their models, theo¬ 
rists find the average over the ensemble of possible outcomes provided by the 
model. 

In inflationary models, as in all passive models, only the initial conditions for 
the perturbations are random. The subsequent evolution is the same for all 
members of the ensemble. For wavelengths higher than the Hubble radius, the 
linear evolution equations for the Fourier components of such perturbations 
have a growing and a decaying solution. The modes corresponding to smaller 
wavelengths have only oscillating solutions. As a consequence, prior to enter¬ 
ing the horizon, each mode undergoes a period of phase “squeezing” which 
leaves it in a highly coherent state by the time it starts to oscillate. Coher¬ 
ence here means that all members of the ensemble, corresponding to the same 
Fourier mode, have the same temporal phase. So even though there is random¬ 
ness involved, as one has to draw random amplitudes for the oscillations of a 
given mode, the time behavior of different members of the ensemble is highly 
correlated. The total power spectrum is an ensemble averaged superposition 
of all Fourier modes and the coherence results in an interference pattern seen 
as the acoustic peaks. 

In contrast, the evolution of the string network is highly non-linear. Cosmic 
strings are expected to move at relativistic speeds, self-intersect and reconnect 
in a chaotic fashion. A consequence of this behavior is that the unequal time 
correlators of the string energy-momentum vanish for time differences larger 
than a certain coherence time. Members of the ensemble corresponding to 
a given mode of perturbations will have random temporal phases with the 
“dice” thrown on average once in each coherence time. The coherence time of 
a realistic string network is rather short. As a result, the interference pattern 

The statistical significance of the observed secondary peaks in the CMB angular 
spectrum is still an open issue, see e.g. Ref. ©. 
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100 1000 100 1000 100 1000 

x c =2.7 x c =6 x c =9 


Fig. 2. Predictions of the toy model of Magueijo et al. (HD for different values of 
parameters x c , the coherence length, and r c , the coherence time; x c oc r}/X c (rj), 
where rj is the conformal time and is the correlation length of the network 

at time ry One can obtain oscillations in the CMB power spectrum by fixing either 
one of the parameters and varying the other. 

in the angular power spectrum is completely washed out. 

Causality manifests itself, first of all, through the initial conditions for the 
string sources, the perturbations in the metric and the densities of different 
particle species. If one assumes that defects are formed by a causal mecha¬ 
nism in an otherwise smooth universe, then the correct initial conditions are 
obtained by setting the components of the stress-energy pseudo-tensor r /iiy to 
zero (Q:0). These are the same as the isocurvature initial conditions (HD- a 
generic prediction of isocurvature models (assuming perfect coherence) is that 
the first acoustic peak is almost completely hidden. The main peak is then 
the second acoustic peak and in flat geometries it appears at l ~ 300 — 400. 
This is due to the fact that after entering the horizon a given Fourier mode of 
the source perturbation requires time to induce perturbations in the photon 
density. 

Causality also implies that no superhorizon correlations in the string energy 
density are allowed. The correlation length of a “realistic” string network is 
normally between 0.1 and 0.4 of the horizon size. 
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An interesting study was performed by Magueijo et al. (]T5|) where they have 
constructed a toy model of defects with two parameters: the coherence length 
and the coherence time. The coherence length was taken to be the scale at 
which the energy density power spectrum of the strings turns from a power 
law decay for large values of k into a white noise at low k. This is essentially 
the scale corresponding to the correlation length of the string network. The 
coherence time was defined in the sense described in the beginning of this 
section, in particular, as the time difference needed for the unequal time- 
correlators to vanish. Their study showed (see Fig. 2) that by accepting any 
value for one of the parameters and varying the other (within the constraints 
imposed by causality) one could reproduce the oscillations in the CMB power 
spectrum. Unfortunately for cosmic strings, at least as we know them today, 
they fall into the parameter range corresponding to the upper right corner in 
Fig. 2. In order to fit the observations, cosmic strings must either be more 
coherent or they have to be stretched over larger distances, which is another 
way of making them more coherent. To understand this, imagine that there 
were only one long straight string stretching across the universe and moving 
with some velocity. The evolution of this string would be linear and the induced 
perturbations in the photon density would be coherent. By increasing the 
correlation length of the string network we would move closer to this limiting 
case of just one long straight string and so the coherence would be enhanced. 


The question of whether or not defects can produce a pattern of the CMB 
power spectrum similar to, and including the peaks of, that produced by the 
adiabatic inflationary models was repeatedly addressed in the literature (§ [T^; 
I7| ; |T8|) . In particular, it was shown (0; O) that one can construct a causal 
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model of active seeds which for certain values of parameters can reproduce the 
oscillations in the CMB spectrum. The main problem today is that current 
realistic models of cosmic strings fall out of the parameter range that is needed 
to fit the observations. 


In addition to purely active or passive models, perturbations could be seeded 
by some combination of the two mechanisms (|5j). Figure 3 (borrowed from 
Bouchet et al.) shows two uncorrelated spectra produced by an inflationary 
ACDM model (dot-dashed line) and a global string model (dashed line), both 
normalized on the COBE data, as well as a weighted sum of the two (solid 
line). The best fit to the latest BOOMERanG, MAXIMA and DASI data 
is obtained with a non-negligible string contribution of 18%. At the moment, 
models combining strings with inflation (|j) or models of varying speed of light 


(p~8| ) are the only ones involving topological defects that to some extent can 
match the observations. One possible way to distinguish their predictions from 
those of inflationary models would be to compute a non-Gaussian statistic, 
such as the CMB bispectrum. 
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Fig. 3. £(£ + 1 )Cg versus £ for three different models [the figure is borrowed from 
Bouchet et ah, Ref. @]. The upper (dot-dashed) line represents the prediction of 
a ACDM model, with ns = 1, Oa = 0.7, fl m = 0.3, fib = 0.05 and h = 0.6. The 
lower dashed line is a typical spectrum produced by global strings. Both lines are 
normalized at the COBE scale (crosses). Combining both curves, with their relative 
contribution as a free parameter, produces the solid curve. The string contribution 
is about 18% of the total. 

3 CMB Bispectrum from strings 

3.1 The string model 


Let us now turn to a computation of the bispectrum of CMB temperature 
anisotropies as generated by simulated cosmic strings. To calculate the sources 
of perturbations we use an updated version of the wiggly string model of 
Ref.@ (see also (|I9|)). 

The network is represented by a collection of uncorrelated, straight string seg¬ 
ments moving with random, uncorrelated velocities, and are produced at some 
early epoch. At every subsequent epoch, a certain fraction of the number of 
segments decays in a way that maintains network scaling. The length of each 
segment at any time is taken to be equal to the correlation length of the 
network which, together with the root mean square velocity of segments, are 
computed from the velocity-dependent one-scale model of Martins and Shel- 
lard (^3|; ^2j ). The positions of segments are drawn from a uniform distribution 
in space and their orientations are chosen from a uniform distribution on a 
two-sphere. 

































According to the one-scale model ([201 ; |2lt |22| ; [Tj|), the string network at any¬ 
time can be characterized by a single length scale, the correlation length L, 
defined by p = p/L 2 , where p is the energy density in the string network. 
Expansion of space stretches the strings without altering their mass per unit 
length. If simply stretched in that manner, strings would soon come to dom¬ 
inate the energy density of the universe. However, at the same time with the 
expansion, long strings reconnect and chop off loops which later decay. These 
two competing effects result in an intermediate steady-state (the so-called 
“scaling solution”) in which the characteristic scale of the string network re¬ 
mains constant relative to the horizon size. 

The total energy of the string network in a volume V at any time r is 

N»L = Vp=^r ( 1 ) 


where N — N{r) is the total number of string segments at that time, V = Voa 3 , 
a is the expansion factor and Vq is a constant simulation volume. From Eq.(l) 
it follows that N = V/L 3 = V 0 /£ 3 , where we have defined i = L/a. The 
comoving length £ is approximately proportional to the conformal time r and 
implies that the number of strings N (r) within the simulation volume Vq falls 
as r -3 . To calculate the CMB anisotropy we need to evolve the string network 
over at least four orders of magnitude in cosmic expansion. Hence we would 
have to start with N > 10 12 string segments in order to have one segment 
left at the present time. Keeping track of such a large number of segments is 
numerically impossible. Fortunately, a way around this difficulty was suggested 
in Ref. (|K[). 

The suggestion is to consolidate all string segments that decay at the same 
epoch. The number of segments that decay at the (discretized) conformal time 
Ti is 


N d (Ti) = V[n(Ti_ i) - n(Ti)] 


( 2 ) 


where n{r) = [£{t )]~ 3 is the number density of strings at time r. Then, with 
the formalism we employ, instead of summing over Ylf= i > 10 12 seg¬ 

ments we now sum over only K segments, making K a parameter. Such a 
simplification only works for computation of two- and three-point correlation 
functions under the assumption that string segments are statistically inde¬ 
pendent. We refer the interested reader to our recent article fl24| ), where we 
explain the procedure in more detail. 

As a final comment, let us note that the simulation model in this form does 
not allow computation of CMB sky maps. This is because the method of 
Ending the two- and three-point functions we describe involves “consolidated” 
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quantities which do not correspond to the energy-moment urn tensor of a real 
string network. These quantities are auxiliary and specially prepared to give 
the correct two- or three-point functions after ensemble averaging. Having 
clarified this, the stress-energy from the simulated string model can be found 
and used to compute the appropriate density perturbations which lead to the 
expected CMB temperature anisotropies. 


3.2 Computation of bispectrum from string realizations 


It is conventional to expand CMB temperature fluctuations into spherical 
harmonics: 


A T/T(e) = J2 ai m Yi m (e). (3) 

lm 


The coefficients ai m can be decomposed into Fourier modes, 

a lm = Cfp(-i) , ‘infd 3 kA l (k) Y ; m (k). (4) 


Given the sources @ M „(k, r) produced by the string model, A;(k) are found by 
solving linearized Einstein-Boltzmann equations and integrating along the line 
of sight (p5|). This procedure can be written as the action of a linear operator 
B[ IV [k) on the source energy-momentum tensor, A/(k) = B[ LV {k)Q r). 

The third moment of aim s can be expressed as 


(a hmi a l2m2 a l3m3 ) = {-i ) h+h+l3 

x Y* m fk 2 )Y* m3 (k 3 ) (A^k,) A i2 (k 2 ) A/ 3 (k 3 )). (5) 

Because of homogeneity we have 

(Aq (k x ) A, 2 (k 2 ) A, 3 (k 3 )) = 5 (3 > (ki + k 2 + k 3 ) P hl2h (k 1; k 2 , k 3 ). (6) 

Given scalar values hi, k 2 , fc 3 , there is a unique (up to an overall rotation) 
triplet of directions k, for which the RHS of Eq. (6) does not vanish. The 
quantity Pip 2 i 3 (ki,k 2 ,k 3 ) is invariant under an overall rotation of all three 
vectors k ? ; and therefore may be equivalently represented by a function of 
scalar values hi, k 2 , k 3 . Hence, we get 


(A h (ki) A; 2 (k 2 ) A i3 (k 3 )) = 5^ (ki + k 2 + k 3 ) p lll2l3 (h, k 2 , k 3 ). (7) 
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The triangle constraint ki + k 2 + k 3 = 0 must be satisfied. Then, in terms of 
the simulation volume Vq, we have 


■^1^3 ^2? ^3) 


Vb 


(A/^ki) A; 2 (k 2 ) A; 3 (k 3 )). 


( 8 ) 


Given an arbitrary direction kx and the magnitudes k \, k 2 and k 3 , the direc¬ 
tions k 2 and k 3 are specified up to rotations by the triangle constraint. Both 
sides of Eq. (8) are function of scalar k t only. The expression on the RHS 
will be computed numerically by averaging over different realizations of the 
string network and over all directions k|. Because of isotropy and since the 
allowed sets of directions k, are planar, it is enough to restrict the numeri¬ 
cal calculation to directions k, in a fixed plane. This significantly reduces the 
computational expense. 

Substituting Eqs. (6) to (8) into (5), Fourier transforming the Dirac delta and 
using the Rayleigh identity, we can perform all angular integrations analyti¬ 
cally and obtain a compact form for the third moment 

(a| imi fll 2 m 2 %m 3 ) = ft™”" 3 / r 2 drb hhh( r )i ( 9 ) 


where, noting ^ fcj 


a Wigner 3j-symbol, we have 


n,m 1 m 2 m 3 _ / ( 2 ^1 + 1 ) ( 2 ^2 + 1 ) ( 2 ^3 + 1 ) 

n hhh - V 47r 



/ 


h l 1 ?. 1 3 


ym 1 m 2 m 3 

and where we have defined, in terms of spherical Bessel functions ji, 


, ( 10 ) 


8 Vq 


b hhh( r ) = ^7J k 1 dkik 2 dk 2 k 3 dk 3 

x 3 h{kir)ji 2 {k 2 r)j h {k 3 r)P h i 2h {k^ k 2 , k 3 ) 


( 11 ) 


Our numerical procedure consists of computing the RHS of Eq. (9) by evaluat¬ 
ing the necessary integrals and averaging over many realizations of the strings 
sources. Readers wishing more details can consult Ref. (j24|). 


From Eq. (9) we easily derive the CMB bispectrum Ci x i 2 i 3 , defined as (|26| ) 




h h h 

yirii m 2 m 3 


Ci 


hhh • 


( 12 ) 
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We now restrict our calculations to the angular bispectrum in the ‘di¬ 

agonal’ case, i.e. h = h = h = l- This is a representative case and, in fact, 
the one most frequently considered in the literature. The power spectrum is 
usually plotted in terms of l(l +1 )Ci which, apart from constant factors, is the 
contribution to the mean squared anisotropy of temperature fluctuations per 
unit logarithmic interval of l. In full analogy with this, the relevant quantity 
to work with in the case of the bispectrum is 


G ul = 1(21 + l) 3 / 2 


f l l l\ 
v o o oy 


Cm. 


(13) 


For large values of the multipole index l, Gm oc l 3 ^ 2 Cm. Note also what happens 
with the 3j-symbols appearing in the definition of the coefficients 7Y)" / 1 J( 2m3 : 
the symbol is absent from the definition of C/p 2 j 3 , while in Eq. (13) 

the symbol (goo) is squared. Hence, there are no remnant oscillations due to 
the alternating sign of (qoo)- 

However, even more important than the value of Cm itself is the relation 
between the bispectrum and the cosmic variance associated with it. In fact, it 
is their comparison that tells us about the observability ‘in principle’ of the 
non-Gaussian signal. The cosmic variance constitutes a theoretical uncertainty 
for all observable quantities and is due to having just one realization of the 
stochastic process, in our case, the CMB sky. 

The way to proceed is to employ an estimator Cq/ 2 z 3 for the bispectrum and 
compute the variance from it. By choosing an unbiased estimator we ensure 
it satisfies C^i^ = (Cw 2 ; 3 ). However, this condition does not isolate a unique 
estimator. The proper way to select the best unbiased estimator is to compute 
the variances of all candidates and choose the one with the smallest value. 
The estimator with this property was computed in Ref. (^) and is 


Ci 


lihh 'y ' 

m\ ,7722,7723 


h I 3 

y rri 1 m 2 U73 


Oil mi Q'l 2 m 2 ^37713 • 


(14) 


The variance of this estimator, assuming a mildly non-Gaussian distribution, 
can be expressed in terms of the angular power spectrum Ci as follows (j2bp 

<7a Cq Ci 2 Ci 3 (1 + ^ij 2 + Si 2 i 3 + <5j 3 ji+ 2<5y 2 Si 2 i 3 ). (15) 
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1 10 100 1000 

multipole Z 

1 /3 

Fig. 4. The CMB angular bispectrum in the ‘diagonal’ case (G ul ) from wiggly cos¬ 
mic strings in a spatially flat model with cosmological parameters = 0.3, 

^baryon = 0.05, = 0.65, and Hubble constant H = 0.65kms _1 Mpc -1 [up¬ 

per panel]. In the lower panel we show the ratio of the signal to theoretical noise 
\Cin/ a Cm | i//;i f° r different multipole indices. Normalization follows from fitting the 
power spectrum to the BOOMERanG and MAXIMA data. 

The theoretical signal-to-noise ratio for the bispectrum (in the diagonal case 
h = h = h = 0 is then given by 

( S/N) h i 2 i 3 = \C h i 2 i 3 /ac w J. (16) 


Incorporating all the specifics of the particular experiment, such as sky cover¬ 
age, angular resolution, etc., will allow us to give an estimate of the particular 
non-Gaussian signature associated with a given active source and, if observ¬ 
able, indicate the appropriate range of multipole /’s where it is best to look 
for it. 


In Fig. 4 we show the results for G)^ [cf. Eq. (13)]. It was calculated using 
the string model with 800 consolidated segments in a flat universe with cold 
dark matter and a cosmological constant. Only the scalar contribution to the 
anisotropy has been included. Vector and tensor contributions are known to 
be relatively insignificant for local cosmic strings and can safely be ignored in 
this model (HI SB The plots are produced using a single realization of the 
string network by averaging over 720 directions of k ? . The comparison of G\(^ 
(or equivalently C)j{’) with its cosmic variance [cf. Eq. (15)] clearly shows that 
the bispectrum (as computed from our cosmic string model) lies hidden in the 
theoretical noise and is therefore undetectable for any given value of l. 


1 The contribution of vector and tensor modes is large in the case of global strings 

d@)- 
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Let us note, however, that in its present stage our string code describes Brow¬ 
nian, wiggly long strings in spite of the fact that long strings are very likely not 
Brownian on the smallest scales. In addition, the presence of small string loops 
(|27|) and gravitational radiation into which they decay were not yet included 
in our model. These are important effects that could, in principle, change our 
predictions for the string-generated CMB bispectrum on very small angular 
scales. 


4 Outlook 


The imprint of cosmic strings on the CMB is a combination of different ef¬ 
fects. Prior to the time of recombination strings induce density and velocity 
fluctuations on the surrounding matter. During the period of last scattering 
these fluctuations are imprinted on the CMB through the Sachs-Wolfe ef¬ 
fect: namely, temperature fluctuations arise because relic photons encounter 
a gravitational potential with spatially dependent depth. In addition to the 
Sachs-Wolfe effect, moving long strings drag the surrounding plasma and pro¬ 
duce velocity fields that cause temperature anisotropies due to Doppler shifts. 
While a string segment by itself is a highly non-Gaussian object, fluctuations 
induced by string segments before recombination are a superposition of effects 
of many random strings stirring the primordial plasma. These fluctuations are 
thus expected to be Gaussian as a result of the central limit theorem. 

As the universe becomes transparent, strings continue to leave their imprint 
on the CMB mainly due to the Kaiser-Stebbins effect. This effect creates line 
discontinuities in the temperature held of photons passing on opposite sides of 
a moving long string. However, this effect can produce non-Gaussian pertur¬ 
bations only on sufficiently small scales. This is because on scales larger than 
the characteristic inter-string separation at the time of the radiation-matter 
equality, the CMB temperature perturbations result from superposition of 
effects of many strings and are likely to be Gaussian. 

Currently, strongest constraints on cosmic string models come from the CMB 
power spectrum data. Predictions of strings-only models are inconsistent with 
data from BOOMERanG, MAXIMA and DASI. However, the data can still be 
fitted reasonably well by “hybrid” models that combine strings and adiabatic 
inflationary perturbations as sources for primordial fluctuations. Presence of 
cosmic strings could, in principle, be distinguished by non-Gaussian statistics. 
The new high-resolution CMB measurements from MAP and PLANCK will 
allow a precise estimate of non-Gaussian signals in the CMB. 

The most straightforward non-Gaussian discriminator is the three-point func¬ 
tion of temperature, or equivalently (in Fourier space) the angular bispectrum. 
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In Ref. (^4|), we have developed a numerical method for computing the CMB 
angular bispectrum from first principles in any active model of structure for¬ 
mation, such as cosmic defects, where the energy-momentum tensor is known 
or can be simulated. The method does not use CMB sky maps and requires a 
moderate amount of computations. We applied this method to the computa¬ 
tion of some relevant components of the bispectrum produced from a model of 
cosmic strings and did not find a non-Gaussian signal that would be observable 
in the forthcoming satellite-based CMB missions. This indicates that either 
the model we used to simulate cosmic strings is inadequate for extracting the 
non-Gaussian signature, or that the angular bispectrum is not an appropriate 
statistic for detecting cosmic strings. 

Avelino et al. (|25| ) applied several non-Gaussian tests to the perturbations 
seeded by cosmic strings. They found the density field distribution to be close 
to Gaussian on scales larger than 1.5(G^/h 2 ) -1 Mpc, where is the fraction 
of cosmological matter density in baryons and CDM combined. Scales this 
small correspond to the multipole index of order l ~ 10 4 . We have not at¬ 
tempted a calculation the CMB bispectrum on these scales because the linear 
approximation is almost guaranteed to fail at such small scales, and because 
of increased computational cost for higher l multipoles. 

The interest in cosmic strings has considerably declined since they have been 
ruled out as the primary seed for structure formation in the universe. This, 
however, is not at all sufficient to rule out strings as an interesting topic of 
research. Studying formation, evolution, properties and cosmological implica¬ 
tions of cosmics strings and other topological defects remains to be both ex¬ 
citing and important. Certain cosmological puzzles, such as baryogenesis and 
ultra-high energy cosmic rays, could be resolved with the help of topological 
defects ([]]). Recent experimental and theoretical progress in understanding 
properties of topological defects in low-temperature condensed matter sys¬ 
tems, such as superfluid helium and superconductors, has led to an intriguing 
possibility of studying Cosmology in the Lab (^). Any new experimental 
constraint on the population and properties of cosmic defects is potentially 
relevant to building viable cosmological and particle physics models. Since 
current and future CMB experiments will provide the most precise cosmolog¬ 
ical data, it is absolutely essential to use them maximally in learning as much 
as we can about cosmic defects. 
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me voy por un camino de labores fecundas, 
a mi muerte yunques taned, callad campanas 0 


5 I’m leaving down a road of fruitful labours, 
anvils toll at my death, hush, bells 
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